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Abstract 

We extend earlier investigations of charge symmetry violation in the valence quark 
distributions of the nucleon, and make similar estimates for the pion. The sensitivity 
of pion-induced Drell-Yan measurements to such effects is then examined. It is shown 
that combinations of tt + and 7r _ data on deuterium and hydrogen are sensitive to these 
violations, and that the pion and nucleon charge symmetry violating terms separate as 
a function of x n and respectively. We estimate the background terms which must be 
evaluated to extract charge symmetry violation. 
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At the present time the flavor structure of the nucleon is a topic of intense interest [1-8]. 
This is largely a consequence of unexpected experimental results, such as the discovery by 
the New Muon Collaboration (NMC) of a violation of the Gottfried sum-rule |§ and the 
so-called "proton spin crisis" [ffOl [11], Q of EMC. There have also been recent theoretical 
calculations of the violation of charge symmetry in the valence quark distributions of the 



nucleon [13, 1 



In most nuclear systems, charge symmetry is obeyed to within about one percent 
15| , so one would expect small charge symmetry violation [CSV] in parton distributions. 



The theoretical calculations suggest that there is a CSV part of the "minority" valence 
quark distributions (d p or u n ), with a slightly smaller violation in the "majority" valence 
distributions (u p or d n ). Although both CSV contributions are rather small in abso- 
lute magnitude, the fractional charge symmetry violation in the minority valence quark 
distributions r min (x) = 2(d p (x) — u n {%)) / \d p \x) + u n (x)) can be large, because at large 
momentum fraction x, d p (x)/u p (x) << 1. Rodionov et al. [14| predicted charge sym- 
metry violation as large as 5 — 10% for the ratio r m i n (x), in the region x > 0.5. The 
relative size of these CSV effects might require a change in the standard notation for 
parton distributions Jl6[ in this region; in addition, Sather [|13| showed that CSV effects 
of this magnitude could significantly alter the value of the Weinberg angle extracted from 
neutral and charged current neutrino interactions. 

Since (in this particular region of Bjorken x) we predict fractional CSV violations as 
large as 5 — 10%, it is important to explore experiments which would be sensitive to the 
relative minority quark distributions in the neutron and proton. Observation of a CSV 
effect at this level would reinforce confidence in our ability to relate quark models to 
measured quark-parton distributions - and hence to use deep inelastic scattering as a real 
probe of the non-perturbative aspects of hadron structure JlT]]. Drell-Yan processes have 
proven to be a particularly useful source of information on the anti-quark distributions in 



nuclei ||18|| . If one uses beams of pions, and concentrates on the region where Bjorken x 
of the target quarks is reasonably large, then the annihilating quarks will predominantly 
come from the nucleon and the antiquarks from the pion. Furthermore, for x > 0.4 
to good approximation the nucleon consists of three valence quarks, and the pion is a 
quark- ant iquark valence pair - in particular, tt + contains a valence d and tt^ a valence 
u. Comparison of Drell-Yan processes induced by tt + and ir~ in this kinematic region 
will provide a good method of separately measuring d and u quark distributions in the 
nucleon. 

Consider the Drell-Yan process in which a quark with momentum fraction x\ in a 
deuteron annihilates with an anti-quark of momentum fraction x<i in a tt + . Provided that 
Xi,X2 > 0.4, to minimize the contribution from sea quarks, this will be the dominant 
process. Neglecting sea quark effects, the Drell-Yan cross section will be proportional to: 

1 
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The corresponding cross section for 7r~D is : 

a^ D ~^(x 1 )+u n (x 1 ))tr(x 2 ), (2) 
so that if we construct the ratio, R^o : 



r°+ D ~-(d p (x 1 )+d n (x l ))d* (x 2 ) 



Arr DY - a DY 
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only charge symmetry violating (CSV) terms contribute. In fact, defining 



Sd = d p -u n , Su = u p -d n , 5T = (f -u« , (4) 

— 7T + 

(and recalling that charge conjugation implies d = d* etc. ) we find that, to first order 
in the small CSV terms, the Drell-Yan ratio becomes: 

) - (£)<*), 

= i^O + iT^), (5) 

Equation (|5|) is quite remarkable in that only CSV quantities enter, and there is a 
separation of the effects associated with the nucleon and the pion. Because is a 

ratio of cross sections one expects a number of systematic errors to disappear - although 
the fact that different beams (tt + and ir~) are involved means that not all such errors 
will cancel. This certainly needs further investigation, since R DY is obtained by almost 
complete cancellation between terms in the numerator. The largest "background" term, 
contributions from nucleon or pion sea, will be estimated later in this letter. However 
we note that Equ. @ is not sensitive to differences between the parton distributions in 
the free nucleon and those in the deuteron |B5L 120, I2T1 E3]. For example, if the parton 



distributions in the deuteron are related to those in the neutron and proton by 

(*) = (1 + (<£(*) + , (6) 

then by inspection equ. (§) will be unchanged. Any correction to the deuteron structure 
functions which affects the proton and neutron terms identically will cancel in R DY . 

It should not be necessary to know absolute fluxes of charged pions to obtain an 
accurate value for R DY . The yield of J/^'s from n + — D and 7r~ — D can be used to 
normalize the relative fluxes, since the J/^'s are predominantly produced via gluon fusion 
processes. The gluon structure functions of the 7r + and 7r~ are identical, so the relative 
yield should be unity to within 1%. 

Next we turn to the predictions for the charge symmetry violating terms, 8d, 5u and 
5d n which appear in equ.(|5]). For the former two there has been an extensive discussion 
by Sather [[0| and Rodionov et al. ||14|| from which there is at least a theoretical consensus 



that the magnitude of 5d is somewhat larger than 5u. This is easy to understand because 
the dominant source of CSV is the mass difference of the residual di-quark pair when one 
quark is hit in the deep-inelastic process. For the minority quark distribution the residual 
di-quark is uu in the proton, and dd in the neutron. Thus, in the difference, d p — u n , the 
up-down mass difference enters twice. Conversely, for the majority quark distributions 
the residual di-quark is a -urf-pair in both proton and neutron, so there is no contribution 
to CSV. 

In Fig. 1(a) we show the predicted CSV terms for the majority and minority quark 
distributions in the nucleon, as a function of x, calculated for the simple MIT bag model 
fl4] , [23| , p4| . There are, of course, more sophisticated quark models available but the 
similarity of the results obtained by Naar and Birse ||25j using the color dielectric model 



suggests that similar results would be obtained in any relativistic model based on confined 
current quarks. The bag model parameters are listed in this Figure. The mean di-quark 
masses are chosen as 600 MeV for the S = case, and 800 MeV for the 5 = 1 case 
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(note that for the minority quark distributions the di-quark is always in an S = 1 state). 
The di-quark mass difference, mad — m uu , is taken to be 4 MeV, a rather well determined 
difference in the bag model. We note that 5u is opposite in sign to Sd and, therefore, 
these two terms add constructively in the Drell-Yan ratio R DY of equ. (|5]). 
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Fig. 1. 

(a) Predicted charge symmetry violation (CSV), calculated using the 
MIT bag model. Dashed curve: "minority" quark CSV term, x5d(x) = 
x (d p (x) — u n (x)); solid curve: "majority" quark CSV term, x5u(x) = 
x (u p (x) — d n (x)). (b) Fractional minority quark CSV term, 5d(x)/d p (x), 
vs. x, as a function of the average di-quark mass = (m uu + m^)/2. 
The di-quark mass difference is fixed at Srrid = vn^ — m uu = 4 MeV. 
From top to bottom, the curves correspond to average di-quark mass 
m d = 850, 830, 810, 790, 770, and 750 MeV. The curves have been 
evolved to Q 2 = 10 GeV 2 . This quantity is the nucleonic CSV term for 
the Drell-Yan ratio of Equ. (0). 
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In Fig. 1(b) we show the fractional change in the minority quark CSV term, 2(d p — 
u n ) I (d p + u n ) vs. x for several values of the intermediate mean di-quark mass. Although 
the precise value of the minority quark CSV changes with mean di-quark mass, the size 
is always roughly the same and the sign is unchanged. This shows that "smearing" the 
mean di-quark mass will not dramatically diminish the magnitude of the minority quark 
CSV term (the mean di-quark mass must be roughly 800 MeV in the 5 = 1 state to give 
the correct N — A mass splitting). 

In Fig. 2 we show the nucleon CSV contribution, R^ D (xi), calculated using the same 
bag model. As is customary in these bag model calculations, this term is calculated at the 
bag model scale (0.5 GeV for R = 0.8 fm) and then evolved to higher Q 2 using the QCD 
evolution equations |[26|| . As the main uncertainty in our calculation is the mean di-quark 
mass (the splitting between S — and S — 1 is kept at 200 MeV [p7|), the results are 
shown for several values of this parameter. In the region 0.4 < x < 0.7, we predict R N 
will be always positive, with a maximum value of about 0.015. For x > 0.7 the struck 
quark has a momentum greater than 1 GeV which is very unlikely in a mean-field model 
like the bag. As a consequence the calculated valence distributions for the bag model 
tends to be significantly smaller than the measured distributions in this region. In these 
circumstances one cannot regard the large, relative charge symmetry violation found in 
this region as being reliable and we prefer not to show it. It would be of particular interest 
to add q — q correlations which are known to play an important role as x — > 1 PHI . 
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Fig. 2. >S 

Contributions from various sources to the nucleonic charge symmetry 

breaking term R^ of equ. (||), evolved to Q 2 = 10 GeV 2 , for pion- 
induced Drell-Yan ratios on deuterons. Unless otherwise specified, pa- 
rameters are: m d = m u = MeV; M (n) = M (p) = 938.27 MeV; 
ft(n) _ jj(p) _ g_g f m . fi md _ mdd _ muu — 4 MeV. Curves represent 

different values of average di-quark mass. Dash-dot: m d = 750 MeV; 
solid: 800 MeV; dash-triple-dot: 850 MeV. 

For the pion, calculations based on the MIT bag model are really not appropriate. In 
particular, the light pion mass means that center of mass corrections are very large, and of 
course bag model calculations do not recognize the pion's Goldstone nature. On the other 
hand the model of Nambu and Jona-Lasinio (NJL) |29| is ideally suited to treating the 
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structure of the pion, and there has been recent work, notably by Toki and collaborators, 



in calculating the structure function of the pion (and other mesons) in this model pO. 51 



The essential element of their calculation was the evaluation of the so-called handbag 
diagram for which the forward Compton amplitude is: 



where 



T, 



fJ,U 



d 4 k 1 



(7) 



T_ = S F (k - q, M^g^T+i^Spik - p - q, M 2 )g 7T g q T^ 5 S F (k - q, Mi), (8) 

represents the contribution with an anti-quark of mass M 2 as spectator to the absorption 
of the photon (of momentum q) by a quark. In equ.(g) g nqq is the pion-quark coupling 
constant, Q the charge of the struck quark, p the momentum of the target meson, and 
Sp the quark propagator. T + is the corresponding term where the quark is a spectator 
and the anti-quark undergoes a hard collision. As in our bag model studies this model 
was used to determine the leading twist structure function at some low scale (0.25 GeV 
in this case), and then evolved to high-Q 2 using the Altarelli-Parisi equations The 
agreement between the existing data and the calculations for the pion and kaon obtained 
in ref.(pOl) was quite impressive. 
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Fig. 3. 

Calculation of the pion CSV contribution R n of equ. (^j), using the model 
of Refs. |30], |3T[ , for average nonstrange quark mass 350 MeV, and mass 
difference — m u = 3 MeV. Dashed curve: R n evaluated at the bag 
scale, fj, = 250 MeV; solid curve: the same quantity evolved to Q 2 = 10 
GeV 2 . 

In the light of the successful application of the NJL model to the structure functions 
of the pion and the kaon, where the dominant parameter is the mass difference between 
the constituent strange and non-strange quarks (assumed to be about 180 MeV), it seems 
natural to use the same model to describe the small difference 5d n (c.f. equ.@) arising 
from the 3 MeV constituent mass difference of the u and dill 511. We have carried out such 
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a calculation. Figure 3 shows the pionic contribution to the CSV Drell-Yan ratio, R n (x 2 ), 
corresponding to this mass difference and an average non-strange quark mass of 350 MeV 
- as used in ref.(|30], |3l|). The dashed curve is the result at the bag scale, and the solid 
curve shows the result evolved to Q 2 = 10 GeV 2 . For x 2 > 0.3, we predict that R* will 
be positive and increase monotonically, reaching a value of about 0.01 at x 2 = 0.8 and 
increasing to about 0.02 as x 2 - ► 1- (We note that there is some ambiguity in translating 
the usual Euclidean cut-off in the NJL model into the cutoff needed for deep inelastic 
scattering. In order to study charge symmetry violation we are particularly concerned to 
start with a model that gives a good description of the normal pion structure function. 
For this reason we have chosen to follow the method used in ref. |3(| rather than ref . [ 31] .) 



Up to this point we have neglected the nucleon and pion sea quark contributions. Since 
the Drell-Yan ratios arising from CSV are very small (viz. Figures 2 and 3), even small 
contributions from sea quarks could make a substantial effect. The dominant contribution 
will arise from interference between one sea quark and one valence quark. Assuming 
charge symmetry for the quark distributions, and using the same form for the quark and 
antiquark distributions in the pion, the sea- valence contribution to the Drell-Yan ratio of 
Eq. (HD has the form 



30(xi,x 2 ) 



a s n Uxi,x 2 ) = ^[2n v (x 2 ) U P(x l ) + 7r s (x 2 )(ul(x 1 ) + dax 1 ))]. (9) 

Unlike the CSV contributions of Eq. (H), the sea- valence term does not separate. In Figure 
4 we show the sea-valence term as a function of Xi and x 2 using recent phenomenolog- 
ical nucleon and pion parton distributions [f32fl . The sea- valence contribution, although 
extremely large at small x, decreases rapidly as x increases. For x > 0.5, the sea-valence 
term is no larger than the CSV "signal". With accurate phenomenological nucleon and 
pion quark distributions, it should be possible to calculate this contribution reasonably 
accurately (the main uncertainty is the magnitude of the pion sea). For smaller values 
(xi ~ x 2 ~ 0.4, where the background dominates, we could use the data to normalize 
the pion sea contribution; we should then be able to predict the sea-valence term rather 
accurately for larger x values, where the CSV contributions become progressively more 
important. We could also exploit the very different dependence on x\ and x 2 of the back- 
ground and CSV terms. We conclude that the CSV terms could be extracted even in 
the presence of a sea- valence "background" . We emphasize that our proposed Drell-Yan 
measurement would constitute the first direct observation of charge symmetry violation 
for these quark distributions. 

The Drell-Yan CSV ratio R®£ of Eq. (|5]) is the sum of the nucleon CSV term of 
Fig. 2 and the pion term of Fig. 3, at the respective values of Bjorken x. Since both 
quantities are positive, they will add to give the experimental ratio. Despite the fact that 
the fractional minority quark CSV term is as large as 10% (c.f. Fig. 1(b)), the nucleonic 
CSV ratio R^ D is more like 1-2%. This is because Sd in equ. (|5|) is divided by u v + d v 
and since d p (x) << u p (x) at large x the nucleon CSV term is significantly diminished. A 
much larger ratio could be obtained by comparing the tt + — p and "7r _ — n" Drell-Yan 
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processes through the ratio : 

A n DY , n DY _ n DY 

^ (Xl ' X2) " Ua°?-v°Y+*°X D )/2- (10) 
To first order in the small CSV quantities, this ratio can be written: 

n ny, , 5d.,f Sd \ . . 

R 7tN {x 1 ,x 2 ) = ^(xi) + I ^+ I (x 2 ), 

= i^(z!) + ir(z 2 ). (ii) 

Once again, the ratio separates completely in xi and and the pion CSV term is identical 
with equ. ([|). However, the nucleon CSV term is much larger - in fact, it is precisely the 
ratio given in Fig. 1(b), so we expect CSV effects at the 5-10 % level for this quantity. 



Fig. 4. 
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CSV terms, but depends on 'EMC 
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EMC term could be significant even for small values of e(x). For x ~ 0.5, where u p /d p ~ 4, 
if e(x) is as large as -0.02 then 5R°^{x = 0.5) ~ 0.10. At larger x the EMC contribution 
could be even bigger, and might conceivably dominate the CSV terms. Since all terms 
(pion and nucleon CSV, and EMC) are predicted to have the same sign in the region 
0.3 < x < 0.8 (we expect e(x) < in this region), the ratio could be as large as 0.3. 
In view of this sensitivity to the EMC term, it is important that accurate calculations be 
carried out of Fermi motion and binding corrections for the deuteron, including possible 
flavor dependence of such corrections. Melnitchouk, Schreiber and Thomas [^, |3^| have 
recently studied the contributions to e(x) in the deuteron. 

We have also calculated the sea- valence contribution to the Drell-Yan ratio R^n ■ 
Relative to the deuteron measurement, we predict a CSV contribution which increases by 
about a factor 5. The sea- valence background also increases by about the same factor. 
So our remarks about the sea-valence background and the CSV "signal" are equally valid 
for these Drell-Yan processes. 

In conclusion, we have shown that by comparing the Drell-Yan yield for tt + and 7r~ on 
nucleons or deuterons, one might be able to extract the charge symmetry violating [CSV] 
parts of both pion and nucleon. We discussed two different linear combinations of 7r + and 
7T~ induced Drell-Yan cross sections, which produce a result directly proportional to the 
CSV terms. Furthermore, we found that the ratio of Drell-Yan cross sections separates 
completely into two terms, one of which (R n (xn)) depends only on the nucleon CSV, 
and the other (R^lx^)) depends on the pion CSV contribution. Thus if this ratio can 
be accurately measured as a function of a; a? and x n , both the nucleon and pion CSV 
terms might be extracted. The largest background should arise from terms involving one 
sea quark and one valence quark. Such contributions, although relatively large, should be 
predicted quite accurately, and may be subtracted off through their very different behavior 
as a function of nucleon and pion x. As the X\ and x-i values of interest for the proposed 
measurements are large (x > 0.5), a beam of 40-50 GeV pions will produce sufficiently 
massive dilepton pairs that the Drell-Yan mechanism is applicable. A flux of more than 
10 9 pions/sec. is desirable, which may mean that the experiment is not feasible until the 
new FNAL Main Ring Injector becomes operable. 

This work was supported by the Australian Research Council. One of the authors 
[JTL] was supported in part by the US NSF under research contract NSF-PHY9 1-08036. 
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